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ABSTRACT 



This paper provides a quantitative account for a recently introduced 
mechanism of mechanical alignment of suprathermally rotating grains. These 
rapidly rotating grains are essentially not susceptible to random torques arising 
from gas-grain collisions, as the time-scales for such torques to have significant 
effect is orders of magnitude greater than the mean time between crossovers. 
Such grains can be aligned either by gaseous torques during short periods of 
crossovers or/and by the difference in the rate at which atoms arrive at grain 
surface arising from difference in orientation of grain in respect to the supersonic 
flow. The latter process, that we call cross-section alignment, is the subject 
of our present paper. We derive expressions for the measure of cross-section 
alignment for oblate grains, and study how this measure depends upon the angle 
between the interstellar magnetic field & the gaseous flow, and upon the grain 
shape. 



Subject headings: magnetic fields - polarization - dust extinction 
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1. Introduction 

In 1949, Hiltner and Hall (see Hall 1949) discovered that the reason for the starlight 
being polarized was orientation of grains of the interstellar dust. More precisely, they 
showed that the polarization was caused by linear dichroism (i.e., by the differential 
extinction of electromagnetic waves of different polarizations), which means that grains in 
the interstellar medium are aligned. 

What makes dust grains (or, at least, a substantial share of them) prefer this or that 
orientation in space? The question is one of considerable nicety and difficulty (see discussion 
in Goodman et al 1995). To search for the possible mechanisms, one should recall several 
basic properties of the interstellar matter (hence, ISM). First, there exist enough evidence 
of paramagnetic nature of a substantial portion of grains (see Whittet 1992). This makes 
them subject to both paramagnetic relaxation 4 and to the Barnett effect, which produces 
magnetic moments of rotating grains. These momenta will then interact with the ubiquitous 
magnetic fields frozen into the interstellar plasma. The interaction will entail precession of 
the moments around magnetic lines. Note that the magnetic moment of a grain is always 
parallel to its vector of angular velocity, which is, roughly speaking, perpendicular to the 
longest dimension of the grain 5 

Strictly speaking, this is valid only for the suprathermal grains, i.e., for those whose 
rotational energy much exceeds that of the Brownian motion. It is these grains that are to 
be of our main concern in this article. We shall address them since a substantial part of 
grains in diffuse clouds should rotate suprathermally (Purcell 1979). 

There may exist several spin-up mechanisms (Purcell 1975, 1979), and the catalytic 



4 Being concerned with mechanical alignment, we do not discuss here paramagnetic 
relaxation, which has its own pros and cons (see Davis & Greenstein 1951, Jones & Spitzer 
1968, Purcell 1969, 1975, 1979, Purcell & Spitzer 1971, Mathis 1986, Roberge, DeGraff & 
Flaherty 1993, Lazarian 1995 a,b,c). 

5 In presence of internal dissipation, a body tends to rotate around the axis of maximal 
moment of inertia. This occurs since the angular momentum is preserved and, for a fixed 
angular momentum, the body will have minimal kinetic energy of rotation, spinning around 
the axis of maximal inertia. As for the inner dissipation, it originates in the following way. 
The angular momentum J of the grain is a product of the inertia tensor by the angular 
velocity co. Generally, co is varying in the coordinate system associated with the grain, and 
so does the magnetic moment \i. As a result, magnetization varies in the grain coordinate 
system and this causes dissipation. 
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one seems to be the most efficient. The essence of this mechanism is as follows. Normally, 
two separate atoms of hydrogen cannot form a molecule because it would have excessive 
energy and, thus, be unstable. The molecule may, however, be successfully produced on 
the surfaces that can trap one hydrogen atom in a potential well 6 ; the second atom will 
react with it, while the excessive energy will be spent for getting the H 2 molecule out of the 
pit, and for acquiring an escape velocity. There may be many such potential wells on the 
surface of each grain (see Tielens & Allamandola 1987, Buch & Zhang 1991). Ejection of the 
nascent molecule will then entail mechanical back reaction, depositing angular momentum 
with the grains and, thus, making them subject to the suprathermal rotation. 

The second fact about the interstellar environment, that we wish to recall, is that the 
ISM is far from being quiescent. On the contrary, various flows are present. Among them 
are stellar winds, collimated outflows from young stellar objects and from accretion discs. 
Moreover, the entire ISM is dominated by supersonic MHD turbulence (see Aron & Max 
1975, Myers 1985, 1987). All these motions of plasma, up to certain reservations, drag the 
dust grains, though with a lag caused by the grain inertia. The lag gives birth to the motion 
of the dust relative to the surrounding gas which, therefore, bombards the grains. Being 
charged grains also spiral about magnetic field (the Larmor frequency is lul) lines and this 
also causes a relative grain-gas motion. A resonance acceleration of grains takes place when 
the Alfvenic frequency becomes equal to As the whole spectrum of frequencies starting 
from a cut-off frequency co c is present the interstellar medium it is possible to show, that 
the acceleration to large relative velocities must be present for sufficiently large grains for 
which lol < <jo c . The relevant detailed study of the dynamics of the charged grains is given 
in Lazarian & Draine (1996, in preparation). 

All in all, if we disregard paramagnetic relaxation, which is marginal for ordinary (i.e., 
without superparamagnetic inclusions) paramagnetic grains in typical ISM conditions, the 
motion of a grain is determined, on the one hand, by its precession around the magnetic 
field line and, on the other hand, by the impacts produced by the flow of the surrounding 
gas. 

The fact that a supersonic flow is capable of providing mechanical alignment was 
discovered by Gold (Gold 1951, 1952). Due to the presence of the precession, the process of 
Gold alignment becomes more complicated, as explained by by Dolginov & Mytrophanov 
(1976), Lazarian (1994) and Roberge, Hanany & Messinger (1995). Both paramagnetic 
and mechanical mechanisms were discussed by the above mentioned authors for thermally- 



6 Such wells (the so-called active sites) may emerge on a grain surface due to, for example, 
inclusions or defects. 



- 5 - 



rotating grains, i.e., for those whose angular momenta are predominantly determined by the 
Brownian motion. The paramagnetic mechanism was generalized by Purcell (1979) for the 
case of suprathermal motions. A similar generalization of the Gold mechanism was carried 
out by Lazarian (1995 d, f). Gold, in his pioneering papers, started out with the obvious 
fact that a flow interacting with a grain produces an angular momentum perpendicular to 
the flow. In case it exceeds the initial angular momentum (associated with the Brownian 
motion), the grain becomes aligned 7 . This consideration is, however, inapplicable to grains 
whose initial angular momenta are suprathermal (i.e. much exceed kT gas ). 

Nonetheless, there exist at least two other, also mechanical, processes that align 
suprathermal dust. The first of them, proposed by Lazarian (1995 d), was based on the 
circumstance that even a suprathermal grain becomes, for a short while, susceptible to the 
gaseous bombardment. This happens at the moment of a crossover, i.e., when the angular 
momentum is minimal. 8 The second process of mechanical suprathermal-grain alignment 
arises from the difference of cross-section of the grain-gas interaction (Lazarian 1995 d,f). 
The essence of this mechanism is as follows. The time between two crossovers depends on 
the rate of gas-grain collisions that is proportional to the gas-grain cross-section. Therefore, 
a grain will tend to spend more time in a position of the minimal cross-section. It is this 
effect that we are going to dwell upon in our article. One may inquire as to which of the 
two mentioned mechanisms will be of leading order in the realistic physical situations. As 
explained in Lazarian (1995d), the answer depends upon the gas velocity. For slow motions 
(< 2 km/s), the random torques caused by hydrogen molecule creation exceed those entailed 
by the gaseous bombardment and, therefore, only marginal alignment may arise during the 
crossovers. Thus, only the second, cross-section mechanism is essential. For high- velocity 
flows, both the crossover and the cross section mechanisms are relevant: their contributions 
to the total effect are of about the same order and, what is especially remarkable, they 
coerce the grain to align in the same direction. It happens as the crossover mechanism 
also tends to minimize the grain cross-section in respect to a gaseous flow (Lazarian 1995 
d). All these estimates are, in fact, valid for oblate grains only. For the prolate ones, the 
cross-section mechanism is weak (Lazarian 1995 d). However, there is evidence (Draine and 
Lee 1984, Hildebrand 1988, Hildebrand & Dragovan 1995) that at least at some regions it 
is oblate grains that dominate polarization. 



7 It is easy to show, that for the mechanical alignment in general and for the Gold 
alignment in particular supersonic motions of grains are necessary. 

8 The crossover happen, for example, due to grain resurfacing (Spitzer & McGlynn 1979) 
caused by accretion: when a monolayer is accreted, active sites emerge on new places. This 
may cause a reversal of the torque direction. 
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Now, we are going to present a detailed calculation of the measure of alignment for 
oblate grains. In the next section, we shall state the topic in the mathematical setting. In 
section 3, the solution to the problem is presented in a general form, while in section 4 some 
particular cases are addressed. The final section, 5, is devoted to a brief discussion of the 
obtained results. 



Our principal goal is to find the measure of grain orientation arising from the 
cross-section mechanism. According to Greenberg (1968), the amount of alignment can be 
expressed by the so-called Rayleigh reduction factor, cr, that relates the statistical properties 
of the grain ensemble to the properties of the observed polarized radiation. 

Let z be the direction of propagation of the observed radiation, and x, y be the 
directions of the maximum and minimum electric wave amplitudes E. Let s be the axis of 
rotation of a spinning grain. Mean absorption cross-sections of (partially aligned, spinning) 
grains for E||cc and E||y will be denoted as C x and C y , respectfully, while those for E||s 
and E_Ls will be named as C(E||s) and C(E_Ls). Then, according to Hildebrand (1988), 
Appendix, 



where 7 is the angle at which the direction, B, of the magnetic field is inclined to the line 
of sight. Then, for spheroidal grains, 



where the angular brackets stand for an ensemble average, and (j> is the angle between 
the magnetic field and the symmetry axis of a grain. As we already mentioned (see the 
comment on page 1), the direction of the magnetic moment will be parallel to the axis of 
symmetry, for oblate spheroidal grains, and will be perpendicular to it for prolate ones. 
So for oblate grains the angle <p will be the one between the magnetic moment and the 
magnetic field. 

The factor a can indeed be considered as a measure of alignment. In absence of 
alignment, ensemble average of cos 2 ip in Eq.(2) become equivalent to averaging over the 
solid angle: (cos 2 ip) = 1/3, and thence a = 0. This gets along with C x being equal to C y . 
For perfect alignment, (cos 2 ip) = 1 and, thus, a is maximal, i.e., unity. 

Thus, our goal now is to calculate a. As explained above, there are indications of the 
polarization effect being mainly produced by oblate grains. So we shall model the grain 



2. Statement of the problem 
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by an oblate spheroid. To compute the average of cos 2 tp over an ensemble, one needs to 
introduce a distribution function p and to execute a normalization 




(3) 



which is to be done in the next section. 



3. 



Calculation of the Rayleigh reduction factor 



For a grain taken at random, its mean time back to crossover, t x , will be of order of the 
mean time between two sequent crossovers, though slightly different. According to Purcell 



provides a good approximation for 0.1 < t^ft^ < 10, being a typical time of existence of 
active sites, and tj, being the rotational damping time (the period during which the gaseous 
friction would stop rotation of a free rotating grain). 

We do not specify here, whether the life-time of active sites is the time necessary for 
accreting a monolayer of refractory material (see Spitzer & McGlynn 1979) or the time of 
poisoning of the active sites (see Lazarian 1995a,f) as our conclusions do not depend on the 
particular process limiting the life of Purcell's rockets. Indeed, what is important for us is 
that, for a chosen grain, is the time of accreting of a particular number of heavy atoms. 
In its turn tj, is proportional to the time of grain colliding with the mass of gas equal to the 
grain mass. 

To understand the essence of the mechanism proposed, consider at first a toy model, 
namely, assume that the axis of suprathermally rotating oblate grain can have only two 
positions, namely perpendicular (_L) and parallel (||) to the axis of a gaseous flux. In this 
model, grain axis stays for the time t x in one position and then undergoes the crossover and 
has equal chances to get either the same or the perpendicular alignment 9 . In the absence of 
the gaseous flux time scales t x \\ and t x ± are the same and, naturally, there is no preferential 
position. However, it is easy to see, that in presence of the gaseous supersonic flux, t x || 
and t x ± are inversely proportional to gas-grain cross sections S n \\ and S n ± respectively. 



9 A complete disorientation is a good model for sufficiently large grains as according to 
eq. (All) in Lazarian (1995a). For small grains the time of complete disorientation can be 
of the order of several t x . This, however, does not change anything in our conclusions. 



(1979) the formula, 




(4) 
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As the time of the crossover is negligible compared with any of t x , we conclude that the 
time averaged probability of finding our grain in positions 1 or || is inversely proportional, 
respectively, to S n ± and S n \\. According to the ergodic hypothesis this probability coincides 
with the ensemble averaged one. 

Returning from the toy model to real grains we observe, that in this case both and 
td are inversely proportional to the rate at which the atoms arrive at the grain surface 
and this rate is proportional to the effective cross-section of the grain, which, in its turn, 
depends upon the angles $, if and $, as on Fig. 1. On this figure, H denotes the magnetic 
field, J stands for the angular momentum, u is the velocity of the gas; angle $, between 
H and, u , varies from to ir; angle if, between H and J, varies from to 7r, while the 
angle ^, between the projections of J and u on the plane orthogonal to H, may vary from 
to 2ir. In the course of precession, the vector J will sweep a cone, of angle if, around 
the magnetic line (Fig. 1). As proven in numerous papers ( Dolginov & Mitrofanov 1976; 
Purcell 1979, Roberge et al. 1993), the period of the precession is many orders less than the 
time t x that has elapsed since the angular velocity of a grain spinning around its own axis 
last crossed through zero. Since the grain preserves the angle if during the time between 
crossovers, and since, by definition, the probability to find a grain at a position given by 
angle if is proportional to the time the grain spends in this position, one can state that this 
probability is proportional to the mean time between two consecutive crossovers first of 
which corresponded to if. 

Damping of spinning is caused both by inelastic collisions and by grain resurfacing that 
changes the positions of active sites. Suppose that on average damping takes place after a 
grain undergoes Ni inelastic collisions, while for a full resurfacing N 2 collisions are sufficient. 
Let N m = max(Ni, N 2 ); then a crossover will happen after approximately N m collisions. 
The frequency of collisions is proportional to the effective cross-section area corresponding 
to a certain position of the grain. Thence, the mean time between 2 consecutive collisions 
is inversely proportional to this effective cross-section. This cross-section varies as a grain 
precesses. Consequently the time t for N m collisions to occur 10 depends on tp: 



where S(t) is cross section of gas-grain interaction, which varies with time. For t much 




(5) 



greater than the period of precession T Eq(5) gives 




(6) 



10 We are happy to acknowledge here, that the elegant derivation of the time between 
crossovers below was pointed to us by Bruce Draine. 
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Inversion of the equation above provides: 

N N 
t = = (7) 

J T nuS(t)dt nuT±S2*S(*)d*' V 1 

where it was taken into account, that the angular velocity of precession lo is equal to 
The probability to find a grain with a given $ is proportional to time t that is stays in this 
position. In fact, Eq.(7) provides the non-normalized distribution function. To normalize 
it, one should integrate (7) also over if, thereby getting for any fixed $ the normalization 
C . Finally, one will get, for the Raleigh reduction factor: 

"=20-2- (8) 

where 



/ dp 

Jo 



sin cp cos 2 — ^ (9) 



and 



C = T dwsmw-z \ (10) 

Jo * ^d* <>(<!>, p,*) 1 ' 

(sin if appears since the integrals in Eqs.(lO) and (9) are to be taken over the solid angle). 

To compute the effective cross-section S'($,(^,^ r ), we approximate the grain by an 
oblate spheroid of a long semi-axis a, and a short semi-axis b (Fig. 2). In Fig. 3, e is the unit 
vector in the direction of the gas flow, while n is the normal to the appropriate cross-section 
of the ellipsoid n (note that this is a regular cross-section, not the effective one which is to 
be perpendicular to the gaseous flow). Let v be the angle between the cross-section and the 

2 2 

axis y. For x = 0, y = pcos v, z = ps'm v, and ^ + = 1, so one can easily express p via 
a, b and v: 

(9 ■ 2 \ — 1/2 

cos z v sin v \ /11N 
— + — j • ("> 

Thus, the area of the figure darkened on Fig. 2 (the regular cross-section) is 

ira ira 
E = Trap = . = . (12) 

V^ + ^ ^cos 2 u + (a/b) 2 sm 2 u 

Now, let a denote the angle between the direction of the gaseous flow and the short axis of 
the spheroid. Then one has, for the unit vector e: e y = sin a, e z = cos a. Similarly, for n: 



11 The circumference of this elliptical cross section corresponds to the points at which the 
gaseous flow is tangential to the grain surface. 
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n y = sin v, n z = cos v. The projection of the area (the regular cross section) to the direction 
will be the effective cross-section we are looking for: 



S = Slnel 



, sin v sin a + cos v cos a 



ira 



cos 2 v + (a/b) 2 sin v 



(13) 



An observer looking at the grain along the vector e will notice an elliptic feature of the 
effective area S. To get rid of the angle v in the latter expression, we would notice that, for 
x = 0, 

y = ay/l - z 2 /b 2 (14) 

and 

tana = — (dy / dz) x=0 = (a/b) 2 tan v (15) 



Plugging this into (13) we obtain an expression for S, in terms of a, b, and a only: 

S = ira 2 cos a ( W tan ' = « a 2 J cos* a(l-(b/a) 2 ) + (b/a) 2 (16) 
<1 + (b/ a) 2 tan 2 a V 



Substitution of this in Eqs. (9) and (10) yields the following expressions for the quantities 
entering formula (8) for the Rayleigh reduction factor: 



I = £ dip sin <p cos 2 <p d*^/cos 2 a(l - ( b/a) 2 ) + (b/a) 2 



C = £ dip sin <p £j d^cos 2 a(l - (b/a) 2 ) + (b/a) 



;i7) 



(18) 



Since the ratio of I and C enter the expression for the Rayleigh reduction factor (see 
Eq. (8)) we omitted in both Eqs. (17) (18) the multipliers l/(7ra 2 ). It follows from Fig.l 
that 

cos a = cos $ cos ip + sin $ sin ip cos ^. (19) 

In section 4 we present analytical results obtained for particular physically important 
cases and numerical computations for other cases. 



4. Results 



4.1. Analytics for particular cases 
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Of practical value is the case of the gaseous flow being parallel to the magnetic 
field. This corresponds to the situation when the dust is driven by radiation pressure (see 
Spitzer 1978). Bearing charge, these grains are coerced to follow magnetic field lines. In 
the reference system associated with the moving dust, the gas will move in an opposite 
direction, i.e., parallel to the magnetic field. (Such motions are typical for bright-star 
vicinities.) In this case, the integrals (9)and (10) can be taken analytically (see Appendix 
A): 



h =0 _ 1 {b/af (a/6) Vl - {b/af - ln(a/6 + ^{a/ bf - 1) 



C$ =0 2 1 -{b/af ln(a/b+ J(a/b) 2 - 1 



The next case of practical importance is that of Alfvenic waves in the interstellar 
plasma. Their interaction with dust is to be discussed elsewhere; here we would just mention 
that according to observations these waves result in supersonic motion of interstellar gas 
(Aron and Max 1975; Myers 1987). The relative direction of grain-gas motion corresponds 
to$ = 7r/2 ( Lazarian 1994). (A similar situation is inherent to the ambipolar diffusion 
mechanism (Roberge & Hanany 1990, Roberge, Hanany & Messinger 1995)). For hypersonic 
flows, it is sufficient to calculate / for exactly $ = 7r/2, and the appropriate asymptotics is 
derived in Appendix B, in assumption of the parameter 

F = l-(b/a) 2 (21) 

being small, i.e., b/a being close to unity. Here follows the asymptotic: 

I\ 112, 

77 «- + — F + F . (22) 

C/j = i,i<i<i 3 45 135 y J 

2 ' 2 a 

We show below, that this approximation is sufficiently accurate for b/a > 2/3. 

Finally, we wish to address the extreme case of dust grains modelled by discs. For an 
absolutely flat disc (b = 0), calculations for $ = ir/2 provide a = 0.25. This is the maximal 
alignment for $ = tt/2. 



4.2. Other cases 



Our analytical studies gave us the grasp of the mechanism under study. To get more 
details we will use direct numerical computation of the integrals involved. For instance, the 
measures of alignment for grains with a number of different axis ratios are shown in Fig. 4 
for $ from to tt/2. Naturally, the difference in cross sections is increasing with the decrease 
of the b/a ratio and therefore the measure of alignment should increase as the result. 
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Fig. 4 proves this tendency. For grains with the axis ratio more than 2/3 the alignment 
is marginal and therefore the influence of the mechanism discussed on the alignment of 
sufficiently spherical grains is negligible. In Lazarian (1994, 1995d) the anisotropy of 
gaseous bombardment was characterized by the parameter s, which is proportional to 
(cos 2 $ — 1/3). It was shown there, that for both the uniform bombardment corresponding 
to (cos 2 $} = 1/3 and the bombardment by an ordered gaseous flux at $ = arccos l/\/3 
produce no alignment of grains. We see that the same result is valid for cross-section 
alignment, namely the alignment vanishes for $ = arccos l/\/3- 

Alignment of grains subjected to Alfvenic waves and ambipolar diffusion is shown in 
Fig. 5. It is easy to see that for grains with the axis ratio greater than 0.1 the measure of 
cross section alignment does not exceed 20%. 

Grains under the radiation pressure follow magnetic field lines and the measure of such 
an alignment is shown in Fig. 6. The efficiency of alignment also decreases with the increase 
of the axis ratio, but for very small axis ratios, i.e. b/a < 0.1 the alignment is more efficient 
as compared with the previous case. 

5. Discussion 

In the article thus far, we have considered alignment of suprathermally rotating oblate 
grains. The alignment is produced in respect to interstellar magnetic field as a grain rapidly 
precesses about this direction. This precession is taking place due to the fact that grains 
are subject to the Barnett effect, i.e., they acquire a magnetic moments as a result of rapid 
spin-up. The alignment itself arises from grain interaction with gaseous flow. 

We have shown that the measure of alignment (the so-called Rayleigh reduction factor) 
is larger for more oblate grains, reaching its maximum for infinitesimally flat disc grains 
(we call them flakes). For example, when $ = (gas flow moves along the magnetic line), 
the said measure 12 reaches —0.5, while for $ = ir/2 (gas flow perpendicular to the field), it 
is 0.25. 

The most striking feature of the discussed mechanism of alignment is that it (unlike 
the Gold alignment) does not depend on the angular momentum that deposited with a 
grain by gaseous flow. It is only the rate, at which atoms arrive at the grain, that matters. 



12 The negative sign of the measure corresponds to preferential alignment of grain long axis 
along to magnetic field, while positive measures correspond to the preferential alignment 
perpendicular to the field. 
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As this rate depends on the grain cross section relative to the gaseous flow, this mechanism 
is called cross-section alignment. To find the direction of the alignment, the rule of thumb 
is that a grain tends to minimize its cross section in respect to the gaseous flow. 

In the vicinity of bright sources grains subjected to the radiation pressure follow 
magnetic field lines. As a result, the cross section mechanism coerce grains to align along 
magnetic field lines. The competing process arises from ubiquitous Alfvenic waves which 
tend to align grains perpendicular to magnetic field lines. In the absence of other alignment 
mechanisms, the interplay between these two effects would determine the polarization 
pattern seen in the ISM. We believe that in reality in different regions of the ISM different 
processes dominate, and one can distinguish these regions only on the basis of comparing 
quantitative theoretical predictions with observations. Such a work for a particular region 
is being done in (Lazarian, Goodman & Myers 1996 in preparation). 

Above we discussed only a very simple model. Indeed, it has been assumed, that the 
flow is hypersonic (so we may disregard stochastic motions of individual atoms), and that 
the direction of the flow does not change in the course of alignment. This is not always 
true. For instance, large-amplitude Alfvenic waves coerce grains to move in such a way 
that the angle between magnetic field and the direction of grain motion is permanently 
changing. In another example, the velocities induced by ambipolar diffusion are usually 
only mildly supersonic and therefore one must take into account the rms velocities of atoms. 
Fortunately, our model can be easily generalized to account for these complications. The 
way for this generalization was discussed in Lazarian (1995d). Without going into details, 
we can exemplify the suggested procedure for the case of rms velocities. In this case it is 
possible to express the angle $ in the following way: 



$ = arctan 

where V{ i = x,y , z are the rms velocity components, while U e are the components of the 
flow velocity. To obtain the corresponding measures of alignment, one needs to integrate 
over the Maxwellian distribution. Obviously enough, only supersonic flows can provide 
substantial alignment; otherwise, if \U\ < \A^, averaging over $ makes the alignment 
negligible. 

In the real world, other alignment processes will also act on the suprathermally rotating 
grains, the crossover (Lazarian 1995d) and paramagnetic (Purcell 1979) mechanisms being 
just two examples. If one knows the measures of alignment for individual processes, i.e. 
crt 1 ) and <j( 2 \ it is possible to find the resulting measure of alignment crS) when the two 



(U X + V X ) 2 + (Uy + Vy) 2 

(u z + v z y 



(23) 
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processes act simultaneously (Lazarian 1995d): 



■(E) 



cr 1 



(i) + ^(1)^(2) + 



(24) 



a 



1 + 2<7( 1 M 2 ) 



Naturally, the above rule can easily be generalized on any number of processes acting 
simultaneously. 

Speaking about the mechanical alignment only, one should bear in mind that, for 
mildly supersonic velocities of grain-gas relative motion, the rms torques associated with H 2 
formation may dominate rms torques arising from atomic impacts during crossovers, i.e., 
"/{{$Jh2) 2 ) > {{S J collision) 2 ) , where 7 _1 is the ratio of the atoms striking the grain to that 
included in the H 2 molecules formed over grain surface. In this case, the alignment measure 
corresponding to crossover mechanism should be diminished by a factor cc/(l + x) where 
x = {(S J collision) 2 ) /["f{{S J H2) 2 )]- If £ <C 1 , the contribution from the crossover alignment 
becomes negligible. Similarly, for ordinary paramagnetic grains in typical ISM conditions, 
the alignment is negligible provided that the spin-up is short-lived (Spitzer & McGlynn 
1979). All in all, there exists a range of conditions where the cross-section alignment is the 
dominant one. 

This work was initiated by Bruce Draine's comments, but would not be possible if not 
for encouragement by Ethan Vishniac. We also indebted to Bruce Draine for pointing to us 
weak points of the original version of the paper and to the referee Wayne Roberge for a 
number of helpful comments. A.L. gratefully acknowledges the support from NASA grants 
NAG5 2773 at the University of Texas at Austin and NAG5 2858 at Princeton University. 
A.L. thanks Alyssa Goodman, Walt Duley, Jim Hough, Phil Myers and Lyman Spitzer for 
illuminating discussions. 



Streaming of grains corresponds to $ = 0. The general expression in this case 
can be simplified and the integration over $ is equivalent to multiplication by 2ir. The 
normalization in this case is equal to 



Appendix 



A. Grains streaming along magnetic field lines 





where 

F = l-(b/a) 2 (A 

and 

E = (b/a) 2 . (A 
Integral (Al) can be easily calculated by substituting z = sin(c^/2) 



irJl - {b/af \ {h/ a) ) 



Similarly the integral / provides: 

W2 cos 2 (,2 sin E ( VF ^ V^+ll 

7o 2tiVF cos 2 tpTTl ~ 2ttF 3 /2 J" 

All in all the ratio I/C is equal to 



_ 1 (6/a) 2 (a/6) yi ~ (6/q) a - Ma/6 + ^{a/bf - 1) 
21-(6/a) 2 ln(a/b+ J(a/b) 2 - 1) 



and this ratio approaches zero for evanescent b/a. 

B. Grain drift perpendicular to magnetic lines 

This case corresponds to $ = tt/2. The integral over $ can be rewritten using the 
notations of the previous section: 

/■tt/2 1 

J = 4/ d^Jsin 2 (p cos 2 (pF + E. (I 
This integral, according to the formula [2.614] in Gradshteyn & Ryzkik (1965), reads: 



J = 4y/E + F sin 2 <pE I , ^ sm y (I 

V^/F sin 2 ip + E ) 

where the complete elliptic integral E of the second kind can be expanded for b/a — > 1 i 
the following way: 



7T 



2 



£(Jfe)«- 1--& 2 -— & 4 - — Jfe 6 -... . (I 



4 64 256 



By means of the latter, the reciprocal former may be presented as 

J" 1 « ^ ( k + -k 3 + -k 5 + —k 7 - 

2vrV^sin^ V 4 64 ^ 256 



where k = \[F sin ip/y F sin 2 cp -\- E. So the normalization will read 



4-7T 



4^ 22 . 19 , 
4: + -F + —F 2 + —F 3 
. 3 30 48 . 



Similarly, for the quantity /: 



7T 



1 2 „ 1 , 155 

- + — F + — F + F 3 

3 15 12 576 



Now we are well equipped to compute the ratio I/C: 
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Fig. 1. — Grain at the position G is subjected to the flux in the direction given by the 
relative velocity u and its short axis is directed along J. The angle ip is measured in the 
plane II which is perpendicular to the direction of magnetic field H. 



Fig. 2. — Cross section of a grain. 
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Fig. 3. — Cross section of a grain. The dotted line corresponds to the crossection S. 
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Fig. 4. — Measure of alignment a as a function of a $, which is the angle between H and u. 
Phi changes from to tt/2. When $ = arccos 1/^3 all the curves pass through zero. The 
solid line corresponds to b/a = 2/3, short dashed line corresponds to b/a = 0.2, while the 
dotted line corresponds to b/a = 0.1. 
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Fig. 5. — Measure of alignment of grains subjected to Alfvenic waves or ambipolar diffusion 
as a function of the axis ratio. This ratio changes from zero (maximal alignment) to one (no 
alignment at all). The dashed line corresponds to the asymptotic obtained for small F. 
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Fig. 6. — Measure of alignment of grains streaming along magnetic field lines as a function 
of the axis ratio. This ratio changes from zero (maximal alignment) to one (no alignment 
at all). The solid line corresponds to the analytical solution found, while the filled hexagons 
corresnpond to numerical computations. 



